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Abstract 



The generalized Heitler-London (GHL) theory provides a straightforward 
way to express the potential energy surface of H3 in terms of Coulomb and 
exchange energies which can be calculated either by perturbation theory or us- 
ing the surface integral method (SIM). By applying the Rayleigh-Schrodinger 
perturbation theory, GHL theory for the quartet spin state of H3 is shown 
to yield results equivalent to the symmetrized Rayleigh-Schrodinger version 
of symmetry adapted perturbation theory (SAPT). This equivalence allows 
a comparison with the corresponding results obtained by the surface integral 
method. The surface integral result calculated with a product of atomic wave 
functions is found to have certain advantages over the perturbation approach. 
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I. INTRODUCTION 



The generalized Heitler-London (GHL) theory provides a useful framework to calculate 
the potential energy surfaces for polyatomic systems Since the potential energy is 

expressed in terms of Coulomb and exchange energies it is possible to systematically separate 
out many-body effects in every single term contributing to the potential energy. In this paper 
some aspects of the three-body exchange effects occurring in H3 are examined in more detail. 

Axilrod, Teller and Muto P were the first to suggest a formula describing the leading 
long range three-body dispersion term for three spherically symmetric atoms. Since then 
the non-additive effects have been intensively studied and several review articles have been 
published @-@|. In the GHL approach the potentials can be decomposed into Coulomb 
and exchange energies, whereas in symmetry adapted perturbation theory (SAPT) these 
interactions are expressed in terms of Coulomb and exchange integrals in the manner first 
introduced by Heitler and London. Recently, SAPT was formulated for the interactions of 
trimers [0] and has been applied to numerical calculations up to third order for the quartet 



spin state of H3 [|1^] and for the helium-trimer []TT| up to third order. Other three-body 
calculations for H3 are based on Heitler-London type calculations |jl2[ and perturbation 
calculations making use of Unsold approximations [0. In the former the splitting into 
Coulomb and exchange part is as pointed out by the author himself not completely rigorous. 

In a previous paper P] analytical results were reported for the doublet as well as for 
the quartet spin state for the H3 system based on the GHL theory. Two kinds of exchange 
energies appear: cyclic exchange energies, where all three electrons are involved, and two- 
body exchange energies in the presence of the respective third atom. The cyclic exchange 
energy of three hydrogen and three helium atoms |n] was calculated using the surface 
integral method (SIM) which was previously applied to two atoms PI, p|, 3, piSHlTl . In a 
forthcoming paper |T8| it will be demonstrated that all exchange energies occurring in the 
H3-system can be calculated either by the surface integral method or by using perturbation 
theory, and the corresponding results for the implicit three-body effect on the two-body 
exchange energies will be derived and compared. 

For H2 it was previously shown that SAPT and GHL are equivalent [|r^ . The purpose of 
this paper is to compare the surface integral method calculations of the three-body effects in 
the exchange energies based on an atomic product wave function with the results of first to 
third order of SAPT which are only available for the quartet spin state of H3 [|1^. In order 
to perform this comparison it is necessary to first prove that the SAPT and GHL theory 
expressions for the energy of the quartet state are equivalent. The results reveal that with 
the zeroth order wave function the surface integral result contains parts of the second order 
SAPT result and is therefore more efficient. 

In Sections || and |T| the basic ideas of the GHL theory and polarization approximation 
are described. In Section the equivalence of the GHL and the symmetrized Rayleigh- 
Schrodinger (SRS) theories is demonstrated order by order. The latter is designated a weak 
symmetry forcing SAPT. Section [V| reviews the surface integral method (SIM). Thereafter 
in Section [V^ the advantages of SIM over the perturbation approach will be demonstrated 
by comparing the numerical results of perturbation theory and SIM. 
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II. GENERALIZED HEITLER LONDON THEORY FOR H3 



The application of generalized Heitler-London theory to H3 was previously discussed in 
Ref. 0. The generalized Heitler-London equation is given by 

HF = Y.^,ng)F (1) 
9 

where F is the localized, i.e. non-symmetrized wave function, T{g) designates a permutation 
operator for the electron coordinates, and eg stands for the Coulomb {g = I) and exchange 
energies {g 7^ I). Applying results from the theory of the symmetric group, the energy 
eigenvalues of the Hamiltonian can be derived. For the Hs-system, the result for the two 
doublet states is 



'^^Eghl = ej - 6123 ±\j\ [(ei2 - £23)^ + (£23 - e,,? + (ei3 - 612)^] (2) 
and for the quartet state 

^^"^Eghl = e/ — ei2 — ^23 ~ ^13 + 2ei23 • (3) 
The remainder of this paper will be concerned only with the quartet state. 

III. POLARIZATION APPROXIMATION AND GENERALIZED 
HEITLER LONDON (GHL) THEORY 

The Born-Oppenheimer non-relativistic Hamiltonian of the three-body system is given 

by 

H = H^ + V (4) 

using 

H' = HI + H% + H'g (5) 
V = Vab + Vbc + Vac (6) 

where H% and Hq are the Hamiltonians of three free hydrogen atoms and Vab? ^bc 
and Vac describe the interaction between atoms A and B and C, as well as A and C, 
respectively. The polarization approximation EO] is based on the equation 



HF = EpF (7) 

where the polarization wave function F and the polarization energy Ep can be written as 
perturbation series 

F = T.^n, (8) 
Ep = Y.^n. (9) 
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The zeroth order polarization wave function 0o is the eigenfunction of the free Hamilto- 
nian and thus is a product of three free hydrogen wave functions. Starting from the 
GHL equation with F chosen as the polarization wave function, Eq. (|1]) together with the 
Hamiltonian Eq. (Q) can be written as 

N 

{H' + V)\T.<Pn) = Y.^gf{g)\Y.ct^n) ■ (10) 

n g n=0 

Forming scalar products with T{g)(f)Q for each group element g 

{f{g)<Po, (H' + V) = ™0o, T.T{g')<Pn) (11) 

n=0 g' n=0 

a system of linear equations can be derived for the Coulomb energy ei as well as for the 
exchange energies eg {g ^ I) in terms of Coulomb integrals J, exchange integrals Kg, and 
overlap integrals Sgi 

Eq + J ^ 6/ + J2g'jtg ^g' S g' -i : g = I 
EoSg +Kg ^ eg + g' ^g g' S g> -1 g ; g ^ I ' 

The following notation for the rath order overlap, Coulomb and exchange integrals was used: 

M 

Sg ■■= E 5." (13) 

n=0 
M 



J ■■=11 J'' (14) 

n=0 

M M 

Kg:=Y.K; = j:Kg, (15) 



n=0 n=l 

where 

S; := {f{g)<Po, <Pn) (16) 

J" := (00, V 0„_i) (17) 

J' = Eo (18) 

K^":=(0o, yf(^7-^)0„_i). (19) 

The equalities Sg-i = Sg and Kg.i = Kg hold. In Ref. |T^ it will be shown how the Coulomb 
and exchange energies can be expressed in terms of Coulomb, exchange and overlap integrals 
and how the order-by-order contributions to the Coulomb and exchange energies can be 
found. 

The convergence properties of the polarization theory have been extensively discussed 
for the case of two hydrogen atoms |21|. For low orders it was shown that the perturbation 



series rapidly converges to the Coulomb energy ||19|, pT|-p3[] though this is not the limit for 
the infinite order expansion. It is assumed that the behavior of this perturbation theory for 
a system of two atoms also roughly holds in the case of three atoms [§,0. Since here we 
are only interested in low orders, especially the first, this expected behavior justifies approx- 
imating the localized wave function via the polarization approximation for three hydrogen 
atoms as well. 
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IV. EQUIVALENCE OF THE GHL AND SRS THEORY FOR QUARTET H3 

In this section the order-by-order equivalence of the complete energy expressions ob- 
tained by using either the GHL or the SRS theory will be demonstrated. Both the GHL and 
SRS theories start with the Hamiltonian Eq. and a zeroth order wave function which is 
a product of three free hydrogen atom wave functions. To demonstrate the equivalence of 
the first order expressions the first order SRS term will be expressed in terms of Coulomb 
and exchange energies. In Eq. (12) of Ref. |^ this term is given by 



3/2 77.1 

^SRS 



<iJo\V {1- t(12) - f (23) - f (13) + f (123) + f (132)) |^o > 



(20) 



which can be expressed with Eqs. (|T6|) to (|19]) as 



3/2 pi 

J^SRS 



where 



-"^23 -"^IS 



-^123 + -^132 



(21) 



qO 
•^12 



•^23 



cO I cO I cO 



^23 



132 



With Eq. (|T2|) it is possible to express the first order contributions as 



-"■123 
-"■132 



e}2 

^23 
^123 



-132 



+ t]S' 
+ e]S\ 
+ e]S' 

^ t/'-'l23 
^ t/'-'l23 




12 

12 


23 

13 



_L.l QO ,1 r-O 
^ t23«J23 -t- fci3^i3 



23^123 



-t- tiQO 



13^123 



^ ^123^123 
^ fcl23*-'23 



-t- tioO 



12^123 



-t- tiqO 



13^123 



-t- t-ioO 



12*-^123 
"12^23 



23*-^123 
t23'-'l3 



1 r-O ,1 r-O 
*^12^13 ^ fc23*JT 



23 '-'12 



ti23*Jl3 
^ ti23*Jl2 

tl3'-'l2 
^ fcl3*-'23 



^ ti320]^23 
^ ^132^13 
' ^132^12 
^ ^132023 

^ ^1320^23 
^ ^123^123 



(22) 

(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



On inserting into Eq. ( plf ) many terms cancel and Eq. (^Tj) is equivalent to the first order 
contribution to Eq. (||) 



3/2 pi 
^SRS 







J' 



123 



^132 



"/ "=12 "=23 '^13 "T "=123 "T "=132 " ^GHL ■ (29) 

The rest of the proof will be done by complete induction. The claim of the induction is the 



equivalence of the GHL and SRS energy expressions up to nth order. From Eq. (12) of |T0 
the general rath-order expression for the interaction energy in SRS theory is found to be 



3/2 7-in 
^SRS 



n-l 



< ^ol (1 - T(12) - f (23) - f (13) + f (123) + f (132)) > 



k=l 



k 

SRS 



< ^ol (1 - T(12) - T(23) - T(13) + T(123) + T(132)) l^^^^^ > 



pol 



No 

n-l 



J -1^12 -"is -"l9S -"■!? 



^■23 



13 



■123 



^32 



E3/2 rrife / nn—k nn—k nn—k i nn—k i Q' 
^SRS \ ^12 ^23 ^13 "T '-'123 "t" : 



n—k\ 
132 ) 



k=l 



(30) 
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where A^'o is given by Eq. 



3/2^; 



GHL 



). Thus it is necessary to prove that 

~ ^12 ti3 -|- ti23 t;^32 



(31) 
(32) 



To perform a proof by induction it is necessary to show that also the (n+ l)st order terms 
of both theories are equal. To do so, the (n+ l)st order of GHL theory is expressed in terms 
of the quantities occurring in SRS theory. This can be achieved by inserting the solutions 
of the set of linear equations Eq. (|T2|) into the complete GHL energy for the Hs-quartet 



state 24 



3/2^ 



GHL 



M 
n=0 



^12 "~ ^23 



GHL 







n=0 



' ^13 + ^123 + ^132 
M 

E - ^12 - e 



23 



1 -s 



12 



s. 



23 



- -I- 1= 
13 ^ ^123 

Sl3 + 'S'i23 + 5*132 



+ e 



132 



(33) 



(34) 



where J, Kg, and Sg have been defined in Eqs. (|T^) to (|T^). To find the expression for the 
(n + l)st order contribution to the energy of the quartet state, the left hand side is first 
multiplied by the denominator 



M 

(E '^'E' 

n=0 



Eq 
M 



M M 

GHL} 1 ~ E('^"2 + 'S'23 + 'S'"3) + E('^123 + 'S'i32y 

n=0 n=0 
M M 

E('^12 + 'S'23 + 'S'^3) + E('^"23 +'^'132) 
n=0 n=0 



+ E - ^r2 - ^23 - ^r3 + ^r23 + 



132J 



n=0 



Collecting terms of (n + l)st order leads to 



(35) 



"'El\\ (1 



cO 



12 



^23 



^13 T '-'I'J?! "T ^1; 



123 



a32j 



jn+1 

+Eo {-S\ 



-"■12 



n+1 
12 



'■23 
ori+l 
'-'23 



on+l 
'-'13 



I on+l I Cfi+l^ 
"T '-'123 "T ^132 ) 



- E "'E^ 



k 

GHL 



_ on+l—k on+l—k on+l— fc 

''^12 ^23 ^13 



ori+l — fe 
^123 



fe=0 

with the result that 



Qn+l—k\ 
•^132 ) 



(36) 



3/2 T^n+l 
GHL 



No 

n 

HE 

k=l 



jn+1 _ ^n+l _ ^n+l _ ^n+l ^ ^n+1 ^ ^n+1 



GHL,k 
3/2 



Q71+1— A; 
''-'12 



'■23 



■13 



^123 



on+l—k 
'-'23 



Q71+1— A; 
'-'13 



Q71+1— A; 
'-'123 



Q71+1— fc^ 
'-'132 J 



(37) 



Using the claim of the proof, which stated that for all orders up to the rath the GHL term 
is equal to the SRS-term, Ey^^''' in the last line can be replaced by '^^'^Eg^j^g^ for all orders 



n. Thus Eq. (|37|) can be transformed into 
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3/2 pn+l 
GH L 



No 



rin+l— fc 
"'-'12 



^23 



^13 



a3 



a32 



(38) 
(39) 



and the equality also holds for the {n + l)st order. Thus the contributions to the energy of 
the Hs-quartet state in the SRS and GHL theories are equal order by order. 



One advantage of the GHL theory is that it permits the calculation of the exchange 
energies by other methods, such as the surface integral method. In Ref. DTOf, the non- 



In Ref. 

additive energy terms of the quartet spin state of H3 have been calculated up to third order. 
The first order terms can be split into a polarization and an exchange part. Since the first 
order polarization energy is pairwise additive, the only non-additive term in first order is 
contained in the exchange term which in Eqs. (23) and (55) of Ref. is given by 

ELh(3, 3) = < ^ol Vab (t(23) + f (13) + f (123) + f (132) - - - 5°23 - ^L) l^o > 



cO 

'-'12 



+ < M Vab (T(12) + T(13) + T(123) + r(132) 
+ < M Vab (t(12) + f (23) + f (123) + f (132) 

which can be expressed in terms of exchange energies as 



cO 
'-'12 



cO 

'-'13 

cO 
•^23 



cO 
'-'123 

cO 
^123 



cO 
^132 

cO 
^132 



1^0 > 
1^0 > 



(40) 



^exch(3,3) 



el23(l 



'S'123) 



6^3(l + ^°0-6t''(l + 5° 



^23 J 



23 J 



e}2(l + 5?2)-4"'(l + 5i°2) 



4(i + ^?3)-4"'(i + 5?3: 



(41) 



This term is also obtained if the pure two-body contributions are subtracted from Eq. (E9| 



V. SURFACE INTEGRAL METHOD (SIM) FOR THE CALCULATION OF 

EXCHANGE ENERGIES 



As shown in Refs. and ||T8| all exchange energies occurring in the GHL-description of 
the H3 system, i.e. the two-body as well as the cyclic exchange energies, can be calculated 
by the surface integral method (SIM). The exchange energy e^g associated with the arbitrary 
group element ^ I is given accordingly by 



-go 



dv 



- {T{go)Ff 



97^1,90 



-1 



dv 



f{go)F] - [f{go)F] V^f} ■ ds 



F{f{gog)F) - {f{go)F){f{g)F) 



(42) 



In order to compare numerical results for three-body exchange effects with the published 
SAPT results for H3 [|lOl, an expression for the non-additive exchange energy has to be ob- 



tained using SIM. The non-additive exchange energy basically contains the cyclic exchange 
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energy and the implicit three-body effects on the two-body exchange energies. As already 
pointed out in Ref. 0] it can be shown that for a choice of the partial volume V such 
that F is localized inside, all quantities occurring in the sum of Eq. ( ^21) go to zero with 
at least a factor of faster than the surface integral itself if all internuclear distances 
are larger or equal to R. This holds for all exchange energies. In a different paper [jT8[ it 
will be shown how to find the implicit three-body effect from the complete surface integral 
expression for the two-body exchange energies. For product wave functions as used here 
the pure two-body part is given by the first line of formula Eq. (^21), i.e. surface integral 
(SI) over denominator. The implicit three-body effect is contained in the second line of Eq. 
(^), i.e. the products of partial overlap integrals with exchange energies. Following the 
same scheme used in the Appendix of Ref. [|14| , these terms can be shown to asymptotically 
go to zero as e~^^ which is faster by a factor of e~^^ than the surface integral (SI) itself. 

Using these results a GHL non-additive exchange energy for the quartet state of H3 
can be defined by simply subtracting the pure two-body contribution from the two-body 
exchange energies in the GHL result for the quartet state Eq. (H) 



(3/2^, 



GHL) exch 



2e 



123 



^12 " ^12 



^23 - ^23 



ei3 



.H2 
= 13 



(43) 



which can be calculated either by SIM or perturbation theory. The first order contribution 
to this non-additive term 



GHL) exch 



2e 



123 



= 12 



H2,l 
= 12 



=23 



H2,l 

^23 



= 13 



= 13 



(44) 



differs from the respective SRS-term Eq. (41) only by overlap integrals that are negligible 
compared to one. 

A comparison of the numerical results of the first order non-additive exchange energy 
Eq. (^) of SRS theory and the GHL term [Eq. (|4^) ] calculated by SIM using the zeroth 



order product wave function F = l/vr^" exp(— — — f3c) is given in Tables | and |I| 
and will be discussed in the next Section. 

In summary, the complete three-body exchange effect in H3, which consists of the cyclic 
exchange energy and the effect of the presence of the third atom on the two-body exchange 
energies, can asymptotically be approximated by the surface integral for the cyclic exchange 
energy. 



VI. RESULTS 

In Tables | and |l| as well as Figures 1 and 2 the numerical results for the first order 
non-additive exchange energy of SRS theory are compared with three different SIM-terms: 
(i) the non-additive exchange energy of GHL theory Eq. (^31), (ii) the cyclic exchange energy 
(complete SIM expression Eq. ( ^2]) with overlaps), (iii) the surface integral (SI) of the cyclic 
exchange energy only (without overlaps). All these quantities have been calculated using the 
zeroth order localized wave function F = l/vr'^/^ exp(— ri^ — — ^30)- Since the exchange 
energies calculated by SIM cannot be given a definite perturbative order (due to the fact that 
only part of the complete space is used in the calculation) the quantity (i) is not expected 
to yield the same numerical results as the first order non-additive exchange energy of SRS 
theory. But since the same zeroth order product wave function was used to calculate both 
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terms it is expected that both quantities exhibit a similar overall behavior in the range of 
parameters studied. 

In Table | results for equilateral triangular geometry of the nuclei ranging between R = A 
and R = 10 atomic units are listed. Generally, all terms calculated by SIM have smaller 
absolute values than the first order perturbative ones. At i? = 4 a.u., the absolute value of 
the complete SIM term Eq. (|^) is 27 % below the SRS result Eq. (|4l|) , the cyclic exchange 
energy is 38 % smaller, and only the surface integral of the cyclic exchange energy is 25 % 
greater in absolute value. At i? = 10 a.u., however, all three quantities calculated by SIM 
are no longer distinguishable and are only 6 % below the SRS result. 

In Table || the results for isosceles triangles with equal sides of length of 6 a.u. and with 
angles 75 varying between 30° and 180° are shown. All quantities except for the surface 
integral without overlaps exhibit a change of sign in the region around 120° and 150°. At 
30°, (i) the absolute value of the SIM term Eq. (^31) is 31 % smaller than the SRS result, 
(ii) the cyclic exchange energy is 41 % smaller, and again (iii) the surface integral of the 
cyclic exchange energy only is 13 % greater in absolute value. At 180° on the other hand, 
only the value for the surface integral has the wrong sign, while both the other terms have 
become indistinguishable and are now 35 % greater in absolute value than the SRS term. 
The differences between the numerical results for the quantities compared in Tables | and 
are, as already pointed out, not due to numerical problems but due to the fact that the 
quantities are different by definition. 

From the Tables it appears that for triangular geometries of the nuclei and internuclear 
distances R > 4 a.u. the first order non-additive exchange energy for the quartet state of 
H3 can be quite well approximated by the surface integral of the cyclic exchange energy. 
This was stated in Ref. [|14[ and has now been explained by the fact that all the SIM 
approximations (see section [V] and in Ref. |T^) hold in this region. 

In Tables |T| and |^ as well as Figures 1 and 2 higher orders of SRS theory are also 
taken into account and compared with the complete GHL non-additive exchange energy Eq. 
(^) in order to show that SIM goes beyond the first order of SRS theory. For equilateral 
triangular geometries of the nuclei and internuclear distances larger than 6 a.u. the results 
of GHL theory lie between the first order SRS term and the sum of the first and second 
order terms, approaching the first order term for increasing distances. At 6 a.u. GHL is 
very close to the first plus second order of SRS, and even at 4 a.u. GHL is only 17 % below 
the total sum up to third order of SRS theory. 

For isosceles structures of the nuclei with equal internuclear distances of 6 a.u. the 
advantage of SIM over the first order SRS theory is even more apparent. Starting at 60°, 
the GHL result is closer to the first plus second order than to the first order SRS term. The 
change of sign occurs for the first order between 120° and 150° whereas for all other terms 
already between 90° and 120°. The differences of the GHL to the first plus second order 
SRS term range from 0.4% at 60° to 33% at 120° and 10% at 180°. At 30° the GHL result 
is again only 16% smaller than the SRS term with the third order term included. 

The advantage of SIM over the perturbative approach is that the surface integral SI is 
easily calculated numerically, and including the partial overlap terms provides part of the 
second order SRS contributions. 
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VII. CONCLUSIONS 



This paper demonstrates how the perturbation series consisting of Coulomb, exchange 
and overlap integrals can be used to express the Coulomb and exchange energies occurring 
in GHL theory. Combining the perturbation series with the GHL theory yields an energy 
expression for the quartet spin state equivalent to that of symmetrized Rayleigh-Schrodinger 



perturbation theory given in ^0 



It is possible to evaluate the exchange energies using the surface integral method (SIM). 
The SIM has the advantage that it derives from a clear physical picture for the exchange 
process in terms of the electrons continuously trading places. For the cyclic exchange energies 
this method has already been described in detail in Ref. [Q, and for the implicit three-body 



effect on the two-body exchange energies it will be shown in Ref. |]TB| . 

The long range behavior of the three-body terms entering the two-body exchange ener- 
gies and of the partial overlap integrals — multiplied by two-body exchange energies in the 
expression for the cyclic exchange energy in Eq. (|^) — indicate that for large internuclear 
separations the surface integral for the cyclic exchange energy is sufficient to describe the 
non-additive contribution to the exchange part of the quartet spin state. The numerical 
results in Tables | and confirm this conclusion. 
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TABLES 



R[ao] 


SRS Eq. @ 


T-il r IT" 1 

EexchlEhl 

GHL Eq. (H) 


2ei23 (SIM) 


2 SI 


4 


-3.83 • 10-3 


-2.79 • 10-3 


-2.39 • 10-3 


-4.21 • 10-3 


5 




-4.31 • 10""^ 


-4.16 • 10-^ 


-5.26 • 10-^ 


6 


-5.90 • 10"^ 


-5.19 • 10-^ 


-5.15 • 10-^ 


-5.70 • 10-5 


7 


-5.88 • 10-6 


-5.32 • 10-6 


-5.31 • 10-6 


-5.55 • 10-6 


8 


-5.33 • 10"^ 


-4.89 • lO-'^ 


-4.89 • lO-'^ 


-4.98 • lO-'^ 


10 


-3.6 • 10-9 


-3.4-10-9 


-3.4 • 10-9 


-3.4 • 10-9 



TABLE I. Comparison of the numerical results for the first order non-additive exchange energy 
of SRS-theory (SRSi Eq. (pl)) with a similar but still different quantity derived from GHL theory 
Eq. (|4^), with the cyclic exchange calculated by SIM (2ei23 (SIM)) including overlaps, and with 
the surface integral SI of the cyclic exchange energy without overlaps (2 SI). The nuclei form 
equilateral triangles with sides of lengths R. 



'jB [degrees] 



SRS Eq. (H) 



ElxchiEh], Rab = Rbc = 6 a.u. 



GHL Eq. (H) 


2ei23 (SIM) 


-2.60 • 10-^ 


-2.23 • lO-'' 


-5.19 • 10-5 


-5.15 • 10-5 


-6.05 • 10-6 


-6.03 • 10-6 


2.61 • 10-^ 


2.60 • 10-"^ 


1.31 • 10-6 


1.30 • 10-6 


1.48 • 10-6 


1.48 • 10-6 



2 SI 



30 
60 
90 
120 
150 
180 



-3.75 • lO-'' 
-5.90 • 10-5 
-7.40 • 10-6 
-3.42 • 10-^ 
8.84 • 10-^ 
1.10 • 10-6 



4.25 


10" 


-4 


5.70 


10" 


-5 


7.95 


10" 


-6 


1.62 


10" 


-6 


5.83 


10" 


-7 


4.10 


10" 


-7 



TABLE II. Comparison of the numerical results of SRS-theory with the same quantities as in 
Table |. The nuclei form isosceles triangles with two sides of lengths Rab = Rbc = ^ a.u., 75 is 
the angle included. 



SRSi Eq. (|ll) 



Eexch [Eh 



SRS2 




SRS3 




-3.60 • 10" 


-3 


-3.34 • 10 


-3 


-5.21 • 10^ 


-5 


-5.03 • 10^ 


-5 


-4.77 ■ 10" 


-6 


-4.62 • 10" 


-6 


-3.71 • 10 


-7 


-3.57-10" 


-7 


-0.7-10- 


9 


-0.7-10- 


9 



GHL Eq. dl) 
-2.79 • 10-3 
-5.19 • 10-5 
-5.32 • 10-6 
-4.89 • lO-'^ 
-3.4 • 10-9 



-3.83 - 10-3 
-5.90 - 10-5 
-5.88 - 10-6 
-5.33 - 10-^ 
-3.6 • 10-9 



TABLE III. Comparison of the numerical results for the non-additive exchange energy in GHL 
theory (GHL Eq. ([4^) with the first order non-additive exchange energy of SRS-theory (SRSi 
Eq. (H)), with the SRS non-additive exchange energy up to second order (SRS2) ||l^ , and with 
up to third order SRS3 |1C] . The nuclei form equilateral triangles with sides of lengths R. 
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7_B [degrees] 



SRSi Eq. (m 



Eexch[Eh], Rab — Rbc — 6 a.u. 
SR,S2 SRS3 



GHL Eq. 



-3.33 


• 10-4 


-3.08 


• 10-4 


-2.60 


• 10-4 


-5.21 


• 10-^ 


-5.03 


• 10-^ 


-5.19 


• 10-5 


-5.67 


• 10-6 


-4.98 


• 10-6 


-6.05 


• 10-6 


3.88 • 


10-7 


9.02 • 


10-7 


2.61 • 


10-7 


1.43 • 


10-6 


1.88 • 


10-6 


1.31 • 


10-6 


1.63 • 


10-6 


2.07 • 


10-6 


1.48 • 


10-6 



30 
60 
90 
120 
150 
180 



-3.75 • 10-4 
-5.90 • 10-5 
-7.40 • 10-6 
-3.42 • 10-7 
8.84 • 10-7 
1.10 • 10-6 



TABLE IV. Comparison of the numerical results of GHL-theory with the same quantities as 

R-BC 



in Table III . The nuclei form isosceles triangles with two sides of lengths Rab 
is the angle included. 



6 a.u., 7b 
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FIGURES 

FIG. 1. Comparison of different orders of the non-additive exchange energy in SRS theory with 
the GHL result (filled triangles) calculated with SIM from Eq. ( |43| ) for equilateral triangles. The 
first order SRS contribution is denoted by circles, and with all terms up to second order by open 
triangles. The stars show twice the surface integral of the cyclic exchange energy. 




4 5 6 7 8 9 10 



internuclear distance R [a.u.] 
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FIG. 2. Comparison of different orders of the non-additive exchange energy in SRS theory 
with the GHL result (filled triangles) calculated with SIM from Eq. ( ^3| ) for isosceles triangles with 
Rab = Rbc = 6 a.u. as a function of the included angle jb- 

The first order SRS contribution is denoted by circles, and with all terms up to second order by 
open triangles. The stars show twice the surface integral of the cyclic exchange energy only. Note 
the change in the energy axis from linear to logarithmic scale. 
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